Finite doping of a one-dimensional charge density wave: solitons vs. Luttinger liquid 
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The effects of doping on a one- dimensional wire in a charge density wave state are studied using 
the density-matrix renormalization group method. We show that for a finite number of extra 
electrons the ground state becomes conducting but the particle density along the wire corresponds 
to a charge density wave with an incommensurate wave number determined by the filling. We find 
that the absence of the translational invariance can be discerned even in the thermodynamic limit, 
as long as the number of doping electrons is finite. Luttinger liquid behavior is reached only for a 
finite change in the electron filling factor, which for an infinite wire corresponds to the addition of an 
infinite number of electrons. In addition to the half filled insulating Mott state and the conducting 
states, we find evidence for subgap states at fillings different from half filling by a single electron or 
hole. Finally, we show that by coupling our system to a quantum dot, one can have a discontinuous 
dependence of its population on the applied gate voltage in the thermodynamic limit, similarly to 
the one predicted for a Luttinger liquid without umklapp processes. 

PACS numbers: 73.21. Hb,71.55.-i,71.45.Lr 
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I. INTRODUCTION 



The transition between a Mott insulator and a metal- 
lic phase in one-dimensional (ID) wires has attracted no- 
table interest in recent yearsi Different types of transi- 
tions are possible, controlled by various physical param- 
eters, such as the electron-electron interaction range and 
its strength. Since the Mott state exists only for commen- 
surate fillings, the electron filling factor inside the wire is 
an important parameter, and it can be varied either by 
controlling the chemical potential, or by doping with a 
finite number of electrons. However, these two methods 
are different: doping with an infinitesimally small num- 
ber of electrons breaks immediately the Mott insulator 
state, while due to the Mott gap, a finite change in the 
chemical potential is required in order to insert the first 
electron and cause the transition 2 . 

The metallic phase resulting from doping is usu- 
ally described by the Tomonaga-Luttinger liquid (TLL) 
theory 2 ^, although for a small finite doping, a large cur- 
vature of the elementary excitation spectrum is expected, 
and deviations from the TLL theory are possible. In the 
TLL, many physical properties are determined by the pa- 
rameter K which describes the interactions. Recently, a 
discontinuity in the occupation of a resonant level which 
is coupled to a TLL with K < 1/2 was predicted by 
Furusaki and Matveev^. However, since the K < 1/2 
regime corresponds to very strong repulsive interactions, 
neither experimental nor numerical evidence for such a 
jump has been obtained so far. 

As a generic model for the wire it is convenient to use 
a tight-binding description of a ID lattice with spinlcss 
electrons. When repulsive interactions between nearest- 
neighbor electrons are considered, the Mott state, taking 
the form of a charge density wave (CDW), occurs for 
strong enough interactions. The Hamiltonian of such a 



system can be written as 



L-l 



H wire = - t y^(ctc 3 - + i + H.c.) 



(1) 
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where I denotes the nearest-neighbor interaction 
strength, and the hopping matrix element between near- 
est neighbors, t, sets the energy scale, ct is the 
creation (annihilation) operator of a spinless electron at 
site j in the L-site wire, and a positive background is 
included in the interaction term. 

The model of Eq. ^ is equivalent to that of the XXZ 
spin 1/2 chain by the Jordan- Wigner transformation 5 . 
The XXZ model is exactly solvable using the Bcthc 
ansatz&i. From this equivalence it is known that for 
a half filled system with periodic boundary conditions, a 
phase transition between a TLL phase and a CDW one 
occurs at / = 2t. Of course, for sufficiently long wires 
the type of the boundary conditions does not change this 
result^. The interaction parameter K is then given by 1 - 



K 



2cos~ 1 (-I/2t)' 



(2) 



so that in the TLL phase, I < 2t and K > 1/2, while for 
/ > 2t, the half filled wire is no longer in the TLL phase, 
but is rather in a CDW state. 

The competition between the TLL and the CDW 
phases is attributed to the presence of umklapp processes 
in this model. For half filling the CDW phase wins over 
the TLL phase once the interactions are strong enough. 
However, when the wire is not exactly half filled, a CDW 
state cannot emerge since it demands a commensurate 
filling, and thus the TLL description is valid even for in- 
teraction values which are greater than 2t. As a result, 
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for strong enough interactions, and sufficiently close to 
half filling, one can then get values of K which are less 
than l/2^i 

Another approach to treat the Hamiltonian of Eq. ([1]) 
is to map it, using bosonization-, into that of the sine- 
Gordon modeli^. In this model, the elementary excita- 
tions are solitons, which carry half an electron charged. 
In the vicinity of the Mott state the value of the TLL 
interaction parameter is K — * j + . The value K = 1/4, 
which cannot be obtained using this model, corresponds 
to the Luther-Emery line, for which the solitons are effec- 
tively non interacting^. For small deviations from half 
filling the solitons are only weakly interacting. 

In this paper we investigate the effects of a finite dop- 
ing on a ID Mott state. We show that although a small 
number of electrons results in an insulator-metal transi- 
tion, the charge distribution in the metallic system with a 
finite number of electrons is not uniform, but rather cor- 
responds to an incommensurate CDW. The added charge 
can also be considered as delocalized solitons in the orig- 
inal commensurate CDW, whose spatial distribution is 
determined by boundary conditions and symmetry con- 
siderations. Nevertheless, the system is conducting, and 
once a bias voltage is applied the solitons are free to move 
and transport charge across the wire. 

An exceptional behavior is exhibited when a lattice 
with an even number of sites L = 2p is occupied by p ± 1 
electrons, i.e., there is a single additional electron (or 
hole) relative to half filling. As will be shown, in this 
case the insulating behavior of the original Mott state is 
retained even though the filling is incommensurate. This 
unusual case will be explained as a result of subgap states 
in the soliton energy spectrum. 

In addition, we show that since the Mott state with a 
small doping allows for the interaction parameter in the 
range 1/4 < K < 1/2, it might exhibit the Furusaki- 
Matveev discontinuity. The appearance of that jump is 
not a foreclosed conclusion, since the model used by Fu- 
rusaki and Matveev does not contain umklapp processes. 
Nevertheless, we find that the jump indeed occurs, and 
show its implications on the solitons present in the wire 
near half filling. 

Since the CDW state is obtained for exactly half fill- 
ing, it is convenient to denote the extra charge by Q = 
n e — L/2, where n e is the number of electrons in the sys- 
tem. As will be shown, the parity of L plays a significant 
role in the behavior when adding the first few electrons. 
Using the finite-size version of the density-matrix renor- 
malization group (DMRG) methodic the Hamiltonian 
H w ire was diagonalized and the ground state was calcu- 
lated for different values of the charge Q and the lattice 
size L. In the current paper we present results obtained 
for / = 3t, which is deep inside the CDW regime (for 
half filling). Nevertheless, other interaction strengths in 
the regime of / > 2t have been checked as well, and were 
found to lead to qualitatively similar results. 

The outline of the paper is as follows: in the next sec- 
tion, Sec. |TT1 we present the charge occupation along 



nearly half filled wires, which demonstrates the presence 
of incommensurate CDW or solitons. In Sec. IIIII we dis- 
cuss the excitation spectrum of the solitons, and the ad- 
dition spectrum of the system. For systems which are ex- 
actly half filled, or doped with only one electron or hole, 
we show that the ground state is insulating; otherwise, it 
is conducting. In Sec. IIVI we discuss the existence of the 
Furusaki-Matveev jump when the doped wire is coupled 
to a resonant level. Finally we conclude in Sec. [V] 



II. CHARGE DISTRIBUTION 

We begin by presenting the occupation of electrons 
along the wire for different numbers of extra electrons Q 
for wires with an odd or an even number of sites (Fig.[T]). 
The Q = case (for an even L) results in a flat distri- 
bution of rij = 1/2 for each lattice site j (dashed line). 
This state is a linear combination of two degenerate CDW 
states, which are eigenstates of the Hamiltonian in the 
thermodynamic limit. These CDW states, which are also 
shown in Fig. [TJ are numerically obtained by applying 
an infinitesimally small potential on the first site, which 
breaks the degeneracy between them£. The population 
of such a CDW state along the wire (near the center) can 
be written as 

rij — n + A cos(27T7y + 0), (3) 

where the amplitude A depends on the interaction 
strength, n — 1/2 is the filling, and is a phase. 

When Q ^ 0, the ground state is no longer degenerate 
and the electron occupation throughout the lead is not 
uniform. Let us start with the Q = 1/2 case (when L is 
odd), which results in a true CDW state, i.e., for each 
even (odd) site the occupation is low (high). Unlike the 
case of the degenerate ground states for Q = 0, the Q = 
1/2 ground state is not coupled by the Hamiltonian to 
the similar CDW ground state of Q = —1/2 due to the 
difference in their total population. 

In order to explore the addition of electrons onto the 
CDW states, one can take the uniform distribution of the 
Q = state as a reference, and investigate the difference 
from it by calculating the accumulated population Arij = 
J2i<j n i ~ i/2< As can be seen in Fig. ® Arij for each 
of the two clean CDW states of Q — (presented in 
solid lines), has an extra charge of e/4 localized at one of 
the wire edges, and a compensating charge (i.e., —e/4) 
localized at the other. The difference between these two 
states is thus a soliton (of charge e/2) located at one 
edge of the wire and an antisoliton (of charge —e/2) at 
the other edge. 

By taking one of the Q = CDW states, and locating 
a soliton at its negatively charged edge, while leaving 
the other edge as it is, one gets a new CDW state in 
which both edges have a positive charge, i.e., a CDW 
state having Q = 1/2. Similarly, placing an antisoliton 
at the edge with the positive charge of the Q = state 
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FIG. 1: (Color online) The electrons occupation in a lattice of 
500 sites (upper panel) and 499 sites (lower panel) for different 
number of electrons. The values near the curves denote the 
extra charge Q which varies between Q = for half filling and 
Q = 3. 



results in the Q = —1/2 CDW state. These states are 
also shown in Fig. (dashed lines). 

For the Q = 1 case (when L is even), adding a local- 
ized soliton near one of the edges is not sufficient, and it 
requires the addition of another charge of e/2. This extra 
charge is obtained by the formation of a delocalized soli- 
ton, centered at the middle of the wire. Further increase 
of the filling above Q = 1/2 (Q = 1) for odd (even) wire 
lengths, results in the addition of two delocalized solitons 
for every electron, so that the total number of delocalized 
solitons is 2Q — 1. 

The delocalized solitons, which carry a charge of e/2 
each, are free to move along the wire. However, when 
no bias voltage is applied, the charge distribution across 
the wire is fixed by the boundaries and by symmetry 
considerations. Practically these constraints lead to a 
density wave with a cosine term similar to that of Eq. ([3]), 
in which the filling factor n is modified. This statement 
is true for much larger values of Q as well. For instance, 
Fig. [3] shows the electron distribution for a 300-site wire 
with Q — 14, i.e., with 27 delocalized solitons. It is easy 
to see that although the amplitude of the oscillations is 
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FIG. 2: (Color online) The extra charge distribution along a 
wire of length L — 300 sites with Q — (solid lines) and of 
L = 299 sites with Q = ±1/2 (dashed lines). 
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FIG. 3: (Color online) The electrons occupation in a lattice of 
300 sites with additional 14 electrons. A fit to Eq. Q (done 
over the middle region) is shown in the inset (symbols). 



reduced, it doesn't vanish and it is rather constant near 
the middle of the wire. Moreover, the oscillations can 
still be fitted to Eq. ([3j using n as a fitting parameter. 
As can be expected, this results in a value of n which 
matches the filling of the wire. 

The non-uniform charge distribution for a fixed num- 
ber of electrons is not a finite-size effect, and survives in 
the thermodynamic limit as well. In Fig. 2] one can see 
a comparison between wires of 300 and 600 sites, with 
Q = 14 and Q — 28. One can see weaker oscillations 
in the Q = 28 case, attributed to a smaller wave length. 
On the other hand, increasing the size of the lead from 
L = 300 to L = 600 while maintaining the number of 
extra electrons Q constant, results in an increase of the 
oscillations. 

One can thus conclude that on the one hand the oscil- 




FIG. 4: (Color online) The electrons occupation near the mid- 
dle of the system in lattices of 300 or 600 sites with additional 
14 or 28 electrons. Results for the 300-site wire are shown as 
a function of 2j. 



FIG. 5: (Color online) The dependence of the addition spec- 
trum on the lattice size L, for different fillings, comparing the 
cases of even and odd sizes and the exact Bethe ansatz results. 
Note the semi-logarithmic scale. 



lations are expected to vanish if a finite filling fraction is 
considered, i.e., in the limit of Q — > oo, L — > oo, keeping 
Q I ' L constant. On the other hand, when a finite number 
of electrons is added, the oscillations are expected to be 
noticed even in the limit of L — > oo. 



III. EXCITATION SPECTRUM 

Up to now we have shown that the density distribution 
of the electrons along the wire, when it is doped by a 
finite number of electrons, is not flat, and preserves some 
features of the Mott state. One can still wonder whether 
these doped states retain the insulating behavior of the 
Mott state or not. To clarify that point we study the size 
dependence of the addition spectrum, defined through 

A 2 (Q) = E (Q + 1) - 2E (Q) + E {Q - 1), (4) 

where Eq(Q) is the ground-state energy of the wire with 
Q electrons above half filling. For a Mott state the limit 
of A2, as L — > 00, should be a finite value, corresponding 
to the gap size, while for a conducting state one expects 
A 2 -> 0. 

The dependence of A 2 (Q) on the wire length L is pre- 
sented in Fig. [5] for even wire sizes (lines) and for odd 
sizes (filled symbols), for different number of electrons 
Q. As can be seen, for Q = and 1/2, the values of 
A 2 converge to the value A 2 (co) (presented by the as- 
terisk symbol at the left margin of the figure), given by 
the Bethe ansatz^. On the other hand, for Q > 3/2 
one gets A 2 — > 0, which indicates that these states are 
conducting. 

A deviation from this intuitive picture appears for 
Q = 1, in which an unexpected gap of A 2 (oo)/2 occurs. 
Nevertheless, this result can be explained by the spec- 
trum of soliton states near half filling. It is known that 
for solitons in the sine-Gordon model with open bound- 
ary conditions there are two degenerate subgap states at 



zero energy, positioned exactly between the conduction 
and the valence bands, each of them localized near one 
of the edges of the system-^. In our case these subgap 
states are thus the localized solitons, which we identified 
as the difference between the two CDW states of Q = 
presented in Fig. [TJ 

Thus the energy spectrum of the solitons is represented 
by valence and conduction bands separated by a gap A, 
while the subgap states exist at A/2 above the valence 
band (see schematic picture in Fig. O . The filling of the 
subgap states with the localized solitons does not change 
the total energy. However, every additional soliton in- 
creases the energy in A/2 (in the limit L — > 00). 
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FIG. 6: (Color online) A schematic picture of the energy 
bands of the solitons according to Ref. [ill The bands are 
separated by energy A, and two discrete states exist at zero 
energy. 

For an odd L, the states with Q = ±1/2, which have 
the same energy due to the particle-hole symmetry of 
the Hamiltonian, differ only in the occupation of the two 
localized solitons. The addition of any other electron is 
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equivalent to the addition of two delocalized solitons, so 
that the energy difference between successive fillings is 
A. Therefore, E {l/2) - E (-l/2) = and E (3/2) - 
E (l/2) = A, resulting in A 2 (l/2) = A. For higher- 
values of Q, one gets Eq(Q + 1) — Eq(Q) = A, so that 
A 2 (Q > 3/2) = 0. These results are summarized in 
Table U 
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TABLE I: Addition spectrum for wires with an odd number 
of sites in the limit L — > oo. 

On the other hand, when L is even, the Q — state 
corresponds to the occupation of only one localized soli- 
ton (or, more precisely, to a linear combination of two 
states, in each of them one edge soliton state is filled and 
the other is empty). The addition of the first extra elec- 
tron fills the additional localized soliton state and a single 
delocalized soliton state, so that Eo(l) — ^o(O) = A/2. 
Every additional electron adds two delocalized solitons, 
thus for Q > 1 one gets E iQ + 1) - E (Q) = A. There- 
fore, as one can see in TableQIl A 2 (0) = A, A 2 (l) = A/2, 
and A 2 (Q > 2) = 0. 
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TABLE II: Addition spectrum for wires with an even number 
of sites in the limit L — > oo. 

Before commencing with the observation of Furusaki- 
Matveev jump we briefly summarize the results so far. 
We have demonstrated the existence of CDW states for 
Q = and Q = ±1/2, and shown that these states are 
insulating. The ground states with Q = ±1 were found 
to be insulating as well, with an excitation gap which is 
half of the Mott gap. States with Q > 1 are conduct- 
ing and thus should be generally described by the TLL 
theory. Nevertheless, we have found that the spatial dis- 
tribution of the electron density is not uniform as in the 
regular TLL picture. Furthermore, this distribution can 
be fitted to an incommensurate CDW form, ft is there- 
fore interesting to explore some other predictions of the 
TLL theory for such states. 

IV. FURUSAKI-MATVEEV JUMP 

The unconventional behavior of the gapless states with 
the non-uniform density can be illustrated by studying 
the discontinuity in the occupation of a resonant level 
coupled to the wirei As mentioned above, once K < 1/2, 



the level is expected to show a jump in its occupation as 
a function of its energy. Nevertheless, one should note 
that umklapp processes, which are an essential ingredi- 
ent in explaining the behavior of our system, were not 
considered in the theoretical framework of Ref. |4j. 

Furthermore, since here the particle distribution in the 
uncoupled wire is not flat, the occurrence of the Furusaki- 
Matveev discontinuity raises an additional question. If 
the electrons occupation profile along the uncoupled wire 
was uniform, one should have observed Friedel oscilla- 
tions in the wire as soon as it is coupled to the resonant 
level. The potential of the resonant level can then be con- 
tinuously changed until the predicted jump in the level 
occupation occurs. At that point, an inversion of the 
Friedel oscillations along the wire is expected. On the 
other hand, once the wire is a Mott state doped by a 
finite number of electrons and the electron density is not 
uniform, the situation is less clear. 

We first calculate the value of K for the doped Mott 
state, and show that it is below 1/2. In order to have 
the same value of K for different lengths of wires, one 
must consider an equal density of additional electrons, 
and we choose to work with doping of Q — L/50, with 
wires of sizes 100 sites and above. K may be obtainec&il 
by the ratio between AE, the energy difference between 
the ground state and the first excited state, and the addi- 
tion spectrum A 2 , since for spinless electrons with open 
boundary conditions AE = ^jf- and A 2 = More ac- 
curate result may be obtained by fitting both A 2 (L) and 
AE(L) to a polynom in l/L, and then obtaining K from 
the ratio of the linear coefficients 18 . Using both methods 
we find that the value of K for our wire is 0.42. 

In order to represent the coupling of a resonant level 
or an impurity of energy eo to the left edge of the wire 
the following term is added to the Hamiltonian: 

H lmv = e a)a-V {a)c 1 + H.c.), (5) 

where a* (a) is a creation (annihilation) operator of an 
electron in the resonant level, and Vq is the hopping ma- 
trix element between the level and the first site of the 
wire. Interaction between the impurity and the wire is 
not considered here since it does not change the result 
qualitatively, but only causes a renormalization of the 
hopping amplitude Vo towards larger values^. 

The level occupation uq as a function of the level en- 
ergy eo is presented in Fig. for different wire lengths 
using Vo — 0.2t. As mentioned above, wires in different 
lengths contain different number of additional electrons 
using Q = L/50, and K = 0.42. When eo is much larger 
than the chemical potential in the wire \x, the impurity is 
almost empty, and the wire contains most of the Q extra 
electrons. On the other hand, for eo <SC /x the impurity 
is almost entirely occupied, and then only Q — 1 extra 
electrons are in the wire. 

Although the population of the impurity is found to 
be continuous for all finite wire lengths studied, it is ex- 
pected to display an abrupt jump in the thermodynamic 
limit, L — > oo. In order to demonstrate this we study the 
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FIG. 7: (Color online) The occupation of an impurity which 
is coupled to one end of ID lattices of different sizes. The 
results shown are for lattice sizes between 100 and 300 (from 
left to right) in steps of 50. In order to compare cases with 
identical values of K, the number of additional electrons is 
taken as Q — L/50, so that K = 0.42. Inset: the slope near 
the point where no = 1/2, which shows a linear dependence 
on the system size. 



dependence of the occupation slope near no = 1/2 on the 
wire length L. As is clearly seen in the inset of Fig. [3 
the slope scales linearly with L, which is a clear sign of 
a first order transition in the thermodynamic limits. 
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FIG. 8: (Color online) The transition of an electron from the 
resonant level into the wire, which leads to the formation of 2 
additional solitons. The results shown are for a level coupled 
to a 150-site wire, with Q — 3. Different curves correspond 
to values of eo between 0.498 (top) to 0.510 (bottom) with 
offsets in the vertical axis for clarity. 

The finite size transition region, in which the electron 



transfers from the resonant level into the wire, is interest- 
ing by itself, since the addition of a single electron to the 
wire is related to the appearance of two additional delo- 
calized solitons in the electron population inside it. As 
can be seen in Fig. as eo increases the electron bound 
to the impurity tunnels into the wire and additional two 
solitons are created in the wire. As mentioned above, 
this is in sharp contrast to the case of a resonant level 
coupled to a TLL having a uniform charge distribution, 
in which only a local change of the Friedel oscillations in 
the charge distribution of the wire is expected. 



V. CONCLUSIONS 

In this paper we have shown that doping a Mott state 
with a finite number of electrons yields states which pre- 
serve the charge modulations of the density wave even 
for an infinite system size. The charge carriers are soli- 
tons, whose spectrum fits the known predictions of the 
sine-Gordon model. 

We have demonstrated how the solitons are added to 
the wire. The first charge of e/2 above half filling results 
in the appearance of a localized soliton at one of the 
wire edges. Each additional charge of e/2 results in an 
additional delocalized soliton, so that for charge Q one 
gets a wire with 2Q — 1 delocalized solitons and a single 
localized one. 

The ground state of the pure Mott state, which is ex- 
actly half filled, is obviously insulating. For wires of odd 
sizes, the ground states for Q = ±1/2 are insulating as 
well, and the excitation spectrum in these three cases 
has the regular Mott gap. Additional insulating states 
with a non-trivial half gap were found for Q = ±1. This 
unusual gap was explained according to the spectrum of 
the soliton states in a ID wire with open boundaries. For 
Q > 1, the excitation spectrum was found to be gapless. 

The unconventional behavior of wires with Q > 1, 
which are gapless but do not preserve translational in- 
variance of charge, is checked against a prediction for a 
TLL wire with an interaction parameter K < 1/2 cou- 
pled to a resonant level. We find that the TLL prediction 
(i.e., the Furusaki-Matveev jump) is indeed obtained also 
for our system, although it involves the creation of two 
additional delocalized solitons in the charge distribution 
of the wire, as opposed to the inversion of Friedel oscil- 
lations in the wire characterizing the conventional TLL 
scenario. As a final remark we note that it may be inter- 
esting to investigate, for this regime of parameters, some 
other TLL predictions as well. 
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